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We investigate a system of harmonically coupled identical nonlinear constituents subject to noise 
in different spatial arrangements. For global coupling we find for infinitely many constituents the 
coexistence of several ergodic components and a bifurcation behaviour like in first order phase tran- 
sitions. These results are compared with simulations for finite systems both for global coupling and 
for nearest neighbour coupling on two- and three-dimensional cubic lattices. The mean-field type 
results for global coupling provide a better understanding of the more complex behaviour in the 
latter case. 



PACS numbers: 02.50.-r, 05.40.+j, 47.20.-k, 47.65. +a 

The influence of noise on nonlinear systems is subject of intense experimental and theoretical investigations jjj. Zero- 
t-H ' dimensional models considering stochastic differential equations for a macroscopic order parameter homogeneous in 
space exhibit noise induced transitions such as transitions between unimodal and bimodal stationary distributions ^,|) , 
cf. also ||] and refs. therein. In real systems like solids and liquid crystals, however, interactions exist between spatially 
distributed constituents. Spatially extended noisy systems described by stochastic partial differential equations are 
difficult to treat analytically, see, e.g. J^L Simulations -though expensive- may provide some guide to a theoretical 
understanding of those systems [|[0J|,| |l0| • 

Models with global coupling of nonlinear noisy constituents are by far easier to investigate and allow even for analytical 
results [p|J^ JTl| . p^ , |l3| , p^| . For example, Shiino JljJ was able to extend the concept of phase transitions to nonequilibrium 

| phenomena described by globally coupled nonlinear oscillators subject to additive noise. More recently, Van den Broeck 
et al. |^J^] demonstrated the appearance of a second order noise induced phase transition in a model with globally 
coupled nonlinear constituents subject to multiplicative and additive noise, which shows no transitions in the absence 
of noise. In this paper we present a model constructed in a spirit similar to jf| which exhibits a first order noise 
induced phase transition connected with a hard onset of the coexisting ergodic components of the system. Varying 

. parameters of the system or of the noise the order of the phase transition may change as we observed previously for 
""ktJ zer0 dimensional models Q. 

Yy- We investigate a system of harmonically coupled identical nonlinear constituents under the influence of noise acting 
simultaneously in additive and multiplicative way. We consider two cases distinguished by the spatial arrangement 
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of the L constituents: The case of global coupling of all components and the case of nearest neighbour coupling on a 
d-dimensional cubic lattice. In the case of global coupling, analytic results are obtained for L — > oo and are compared 
with simulations for L = 100 and L = 1000, respectively. Furthermore, simulations were carried out for the case of 
nearest neighbour coupling in d = 2 (L = 100 x 100) and d = 3 (L = 20 x 20 x 20). In these cases the results for 
^ . global coupling can be considered as mean field approximation. 

The variables Xi of the individual constituents at the lattice sites i obey the following stochastic differential equations 
in the Stratonovich sense 



Xi = f(xi) + g(xi) & - ^ ( Xi ~ x i) ■ W 

3'eAT(i) 

Here Af(i) denotes the set of involved neighbours of site i and N = #J\f(i) is equal to L — 1 in the case of global 
coupling and to 2d in the case of nearest neighbour coupling. The parameter D controls the strength of the spatial 
interactions. The represent zero mean spatially uncorrelated Gaussian white noise at point i with autocorrelation 
function 

<£<(*) &(* , )> = " 2 M(*-*'), (2) 

where a 2 is the noise strength. For nearest neighbour coupling and suitable chosen parameters, Eq. (|l|) can be 
considered as discretized version of a stochastic partial differential equation with diffusive coupling. 
The stationary Fokker-Planck equation for the probability density of a;, reads 

= ljx(~^ Xi ) + ~Fr 13 ( Xi ~ fol 35 *)) + ^dix^^-gixS^Psixi) , (3) 
1 V jeAf(i) 1 J 
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where {xj\x{) = J dxjXjP s (xj\xi) is the steady state conditional average of Xj, j £ A/"(i) given Xi at site i. 
For the case of global coupling, fluctuations disappear in the average 1/(L — 1) Ylj£j\f(i)( x j\ x i) if Z — > oo. Considering 
the class of solutions for which this expression is independent of lattice site i we can replace it by the steady state 
mean value (x) which is self-consistently determined by 

dxxP s (x,{x)) = F({x)) . (4) 

Obviously, for finite lattices or coupling to a finite subset of neighbours this replacement represents a mean field 
approximation. 

Following Shiino Jl|| one obtains the same results by replacing in (|l|) the spatial average l/^EjeA/Yt) x j by tne 
statistical average (x). 

In this paper we consider a simple model, for which a nontrivial solution of (^) is not emerging from zero but appears 
with a jump to a nonzero value at a critical value of the control parameter. The model is specified by 

f(x) = ax + x 3 - x 5 , g(x) = 1 + x 2 . (5) 

For the model with global coupling the stationary probability density is 



J i 

P s (x,(x)) oc (l + x 2 )° 2 exp<jcr" 2 



. , D - a + 2 



x H - h D (x) [ x + arctanfx) 

1 + X z iii.. 



(6) 



Without spatial coupling (D — 0) the model shows the following bifurcation behaviour: 

In the deterministic case (a 2 = 0) the stationary solution undergoes a subcritical bifurcation at a c — 0. The multi- 
plicative noise shifts the bifurcation threshold of the maximum of the stationary probability density to a c = a 2 . For 
weak noise (a 2 < 1) the bifurcation is subcritical whereas for a 2 > 1 it is supercritical. 

In the following we mainly restrict ourselves for the sake of convenience to the case a 2 — 1 where the noise intensity 
is just sufficient to produce a change from the deterministic subcritical bifurcation into a supercritical bifurcation. 
Then the extreme values of the stationary probability density are x s t = and, if a > 1, x s t — ±vo— 1. 
The global coupling (D ^ 0) favours a coherent behaviour of the spatially distributed components which is -in a 
sense- an effect opposite to the noise and will 'restore' the subcritical bifurcation. Hence we expect to observe a first 
order nonequilibrium phase transition. 

The bifurcation behaviour of (x) is governed by the self-consistency condition (jfy. Since for our model F ((x)) is an 
odd function of (x) we always have the solution (x) = 0. Moreover, pairs of new stable and unstable nonzero solutions 
may occur in certain parameter ranges. Only stable solutions can be observed in simulations. We remark without 
proof that F' ((x)) < 1 is sufficient for stability (cf. |[3)). The existence of more than one stable solution leads to the 
existence of several corresponding stationary probability densities P s {x, (x)). Therefore, a phase transition breaking 
the ergodicity of the system is expected in the case of global coupling. 

The typical behaviour of F ((x)) for our model is sketched in Figure 1. Whereas the model investigated by Van den 
Broeck et al. || exhibits a second order phase transition, we find a first order phase transition connected with a hard 
onset of the nontrivial stable solution of the self-consistency condition for our model in a certain parameter range 
where the model without spatial coupling exhibits only the trivial solution (x) =0. 

The phase diagram given in Figure 2 confirms the intuitive picture drawn above: The spatial coupling favours coherent 
behaviour of the components acting thus opposite to the noise. The critical value of a is reduced with increasing 
coupling strength D and above a critical strength of D the transition is of first order. 

Figure 3 shows the different solutions of the self-consistency equation (^) for the order parameter (x) leading to 
different ergodic components of the system as a function of the parameter a, the spatial coupling constant D, and the 
noise strength a 2 , respectively. In all cases a hard onset of the nontrivial stable solutions can be observed. 
It is instructive to compare the results obtained by solution of the self-consistency condition (Q) with simulations 
on finite globally coupled systems of different size. We consider a parameter set (a 2 = 1,D = 25, a = —1.5) where 
three different stable solutions of (Q) exist. For small systems (L = 100) the ergodicity breaking is not perfect. We 
still observe a few transitions between the different 'ergodic components' due to large fluctuations of (x(t)) around 
its stationary values. The trajectory of the spatial average (x(t))i = Y2j—x x i(t) is shown in Figure 4. For 

larger systems (L = 1000) the fluctuations become smaller and the system remains very long inside one of the ergodic 
components. In that case, there are practically no transitions. The initial conditions determine which of the ergodic 
components is selected. 
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We also performed simulations of the stochastic differential equation (|l|) for nearest neighbour coupling on a 3- 
dimensional cubic lattice with L = 20 x 20 x 20 sites and on a 2-dimensional square lattice with L = 100 x 100 sites. 
We used the same set of parameters as in the simulations for global coupling. Qualitatively one gets a behaviour very 
similar to the case of global coupling. In Figure 5 we compare the probability density P s {x) at one lattice site given by 
for the case of global coupling with the results of simulations for nearest neighbour coupling on the 3-dimensional 
lattice. Although for the finite system there is no perfect separation into different ergodic components, the trajectories 
remain very long in the corresponding 'basin of attraction'. The histograms obtained by sampling those trajectories 
follow very closely the probability densities of the ergodic components for global coupling. The value of x s for nearest 
neighbour coupling is about 10 percent smaller than for global coupling. 

Simulations on a 2-dimensional square lattice (L = 100 x 100) exhibit a similar qualitative behaviour, the quantitative 
agreement with the results of the globally coupled model is -as has to be expected- less satisfying. Figure 6 compares 
the order parameter as a function of a and a 2 for the 2-dimensional lattice with the results for global coupling. Al- 
though the bistable region is smaller than in the case of global coupling it remains no doubt that it exists. We remark 
that the fluctuations (indicated by the error bars in Fig. 6) are larger for the states with (x) ^ than those for (x) = 
which is a clear indication of the multiplicative nature of the driving process. As in the case of global coupling, for 
the L = 100 x 100 system the trajectory stays inside the ergodic component selected by the initial condition for a very 
long time. No jumps were observed in our simulations running typically over a time t — 5000. Jumps between the 
ergodic components induced by large fluctuations are observed in simulations with smaller systems. With increasing 
size they become less frequent, cf. Fig. 7. 

In this paper we investigated a model which exhibits a first order noncquilibrium phase transition due to a hard onset 
of the coexistence of several stable ergodic components of the system. Other models which we investigate at present 
exhibit the same behaviour. In a different context, a system of coupled Duffing oscillators was used to caricature a 
liquid to crystal transition Jl6| ]. Mean field theory yields a first order nonequilibrium transition, which is preserved 
including additive noise. We also found changes from second order to first order transitions by tuning parameters of 
the noise or the system. In any case, both the nonlinear terms and the interplay between deterministic and stochastic 
effects determine the order of the transition. In previous work fl|7j-|s|.pHTc|] only second order noise-induced nonequi- 
librium phase transitions have been observed. 

Our results may be of interest in the context of experimental investigations in electrohydrodynamic convection in 
nematic liquid crystals subject to thermal fluctuations (additive noise) and/or an external stochastic voltage (mul- 
tiplicative noise). There are experimental hints jl7| that the first transition from the homogeneous state to the 
structured state might be weakly hysteretic although the deterministic theory predicts a supercritical bifurcation. 
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NOTE ADDED 



After completion of this work we got knowledge of a paper by S. Kim, S.H. Park, and OS. Ryn (Phys. Rev. Lett. 
78, 1616 (1997)) where a first order noise induced transition on a different system of globally coupled oscillators is 
described. 
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FIG. 1. Solution of the self-consistency equation (Q), F((x)) = (x), in three typical cases, (x) = is always solution; 
in the case of the dashed-dotted line it is the only solution. In the case considered by Van den Broeck et al. Q (dashed 
line) we have two stable solutions (x) = ±x s (full circle), (x) — is unstable. In the case considered here (solid line) 
we have besides the stable solution (x) = a pair of unstable solutions (x) — ±x u (empty circle) and a pair of stable 
solutions (x) ~ ±x s (full circle). In contrast to the former case, in the latter case the nontrivial solutions do not 
emerge continuously from (x) — but appear with nonzero value at the critical value of the control parameter. This 
indicates a first order nonequilibrium phase transition. 




FIG. 2. Phase diagram in the case of global coupling for a 2 — 1. For small D we have a second order transition. The 
spatial coupling favours a coherent behaviour of the constituents, acting thus opposite to the noise. With increasing 
coupling strength D the critical value of a is reduced and above a critical strength of D the first order transition 
of the model without noise and spatial coupling is 'restored'. The solid (dashed) lines denote a first (second) order 
nonequilibrium phase transition. The number of ergodic components is three in the shadowed region, two in the 
region above and one in the region below. Hysteresis appears in the shadowed region. 
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FIG. 3. Stable solutions of the self-consistency equation ([|) as function of the control parameters (solid lines) 
determine the ergodic components. Unstable solutions (not shown here) cannot be observed in simulations. The 
maxima of the stationary probability density (dashed line) of the corresponding ergodic components (|^) exhibit 
qualitatively the same behaviour. In all the diagrams we observe at the critical value a hard onset of the nontrivial 
stable solution corresponding to a first order transition. In (c) a reentrant behaviour is found similar to that in ||. 
Parameters are a 1 = 1, D = 25 in (a), a 1 = 1, a = —1.5 in (b), and a = —1.5, D = 25 in (c). 
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FIG. 4. Trajectory of the spatial average (x(t)) l = (l/£) x i(t) f° r the case of global coupling (L — 100, a 2 = 
1, a = —1.48, D = 30). The trajectory fluctuates preferably around the mean values x s — and x s — ±0.94, sometimes 
large fluctuations lead to jumps between the 'ergodic components'. 




FIG. 5. Probability densities for Xi(t) at arbitrary i for the case of global coupling (solid lines) as given by (^). The 
ergodic components correspond to (x) — 0, x s = 0.94 [a 2 = 1, D = 25, a = —1.5). These results are compared with 
simulations for the 3-dimensional cubic lattice (L = 20 x 20 x 20) with nearest neighbour coupling. The probability 
densities (o and o) are obtained by sampling 20000 equidistant points within a trajectory of length t = 10000 near 
(x) = and x s — 0.85 separately. The plot indicates that the globally coupled system gives a good idea of the 
qualitative and quantitative behaviour of the system with nearest neighbour coupling. 
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FIG. 6. Comparison of the order parameter (x) obtained by simulation for a 2-dimcnsional square lattice of size 
L = 100 x 100 with the results for the globally coupled model (thick solid line) for D = 30. The diamonds denote the 
average of Xi(t) over all lattice sites and over a time span of order 100 during which no jumps between the ergodic 

— 1/2 

components occur. The error bars indicate the time average over the standard deviation (l/L ^(ar^i) — (x(t)) L ) 2 ^ 
(a) and (b) show the dependence on the control parameter a (a 2 = 1) and the noise strength a 2 (a = —1.5), 
respectively. The coexistence of the solutions with (x) ^ and (x) = over a range of parameters as for the case of 
global coupling is obvious. 
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FIG. 7. Trajectories of the spatial average (x(t)) l for the 2-dimcnsional square lattice of size (a) L = 10 x 10 and (b) 
L = 18 x 18 for D = 30, cr 2 = 1, a = —1.48. For the smaller system frequent jumps between the 'ergodic components' 
are observed; with increasing size of the system these events are rarefied. Already for a size of L = 100 x 100 no 
jumps were observed in time spans of order 5000. 
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